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1. Introduction 

The purpose of this paper is to explore the geometry of the inextensible string 
(or whip) in Euclidean space , and its application as an alternative geometry in 
shape recognition to the geometries proposed by Michor-Mumford [MM 11 IMM2j , 
Younes et al. |YMSM| . Klassen et al. |KSMJ| . and others. Generally N ^ 2, 
although we will assume whenever convenient that N = 2; this simplifies formulas 
but does not substantially change any of the results. 

In the absence of external forces, the string is a geodesic motion in the space 
of unit-speed curves in R^; hence it is another of the examples of partial dif- 
ferential equations arising as geodesic motion on an infinite-dimensional manifold 
which have been discovered in the wake of Arnold's jX] approach to hydrodynamics 
(other examples include the Korteweg-deVries equation [OKj . the Camassa-Holm 
equation [M|, and the Hunter-Saxton equation [KM]). The equations of motion are 
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with initial conditions 7?(0, s) = 7(5) and 7?t(0, s) = i«(s), assumed to satisfy the 
compatibility conditions |7'(s)| = 1 and (7'(s), w'(s)) = 0. Here t is the time 
parameter, s is the length parameter along the curve, and a is the tension. We 



2 



STEPHEN C. PRESTON 



will suppose the curves have one fixed and one free endpoint. We will discuss other 
boundary conditions in Appendix |Al 

As shown in the author's companion paper jPrj (where a slightly different nota- 
tion was used to make the estimates more convenient), the easiest way to handle the 
fixed point is to extend the curve through the fixed point by oddness to get a curve 
with two free endpoints. So we have a curve 77: [0,T) x [—1, 1] satisfying 
ri{—s) = —ri{s), and the boundary conditions for © are then cr(— 1) = cr(l) = 0. We 
then automatically have, for any sufficiently smooth solution, that (7(—s) = cr(s), 
so that all even derivatives of rj and all odd derivatives of a vanish at s = 0. Using 
a spatial discretization (the method of lines), the author proved local existence and 
uniqueness in the weighted energy norm 

m pi 

(3) Em = J] \ {l-sy\divt\' + {l-sy+'\dt'v\'ds 

for m ^ 3. Our aim in this paper is to explore the geometric interpretation of 
this result, and on the way we will obtain a result on the dependence of solutions 
on the initial conditions. Specifically, we show that for a fixed 7, the solution is 
differcntiable but not as a function of w. 

The system ([T|)-([2|) have been known and studied for hundreds of years, although 
only recently has there been a rigorous proof of well-posedness for the full nonlinear 
system (see [Pr]). Thess et al. jTZNj studied these equations on the circle as 
a toy model of hydro dynamical blowup, and one of our motivations in studying 
the geometrical aspects was to see just how far this analogy goes (comparing to 
Arnold's geometrical approach to hydrodynamics 0). See also Serre [S] and Reeken 
[RH IR21 IRS) for analytical aspects of these equations, and references cited in |Pr| 
for physical discussions of their properties. 

First, in Section [2] we define precisely the manifold structure on the space of 
curves and show that the inextensible curves form a smooth submanifold. Since 
we want to use the implicit function theorem to do this, we want to work on a 
Hilbert manifold, and the result ([3]) suggests the appropriate Sobolcv topology is 
the weighted space 

(4) /C™ = {t]-. [-1, 1]^R^ :Rr] = ij and \\r]\\j^j < 00 for ^ j ^ m} 

where R is the odd reflection operator {Rri){s) = — ?;(— s) (which is bounded in any 
weighted Sobolev topology) and the weighted energy norm is defined by 

(5) Mik = f^{i-sy\d':7j{s)fds. 

We will first prove that the configuration space 

(6) A"' = {?7 6 /C'" : \t]'\ = 1 and i?ry = 77} 

is a smooth submanifold of /C™ for m ^ 4, with tangent space given by 

(7) T^A"" = {f 6 /C™ : (v', -f'y = and Rv = v}. 

(The space A"^ is a smooth manifold as long as m ^ 2, but it fails to be a subman- 
ifold of /C™ if m = 3 or m = 4.) 
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In Section [3] we define a weak Riemannian metric on /C™ (and hence A'"'), given 
for vector fields u and v along a curve rj, by 



Although we have a smooth metric on a smooth manifold, the Lcvi-Civita connec- 
tion is not smooth, and thus the geodesic equation is not an ordinary differential 
equation. (This is typical behavior for a weak metric on an infinite-dimensional 
manifold; the Levi-Civita connection is unique if it exists, but it is not even guar- 
anteed to exist if the Riemannian metric does not generate the topology of the 
manifold; see |EM| .) This is reflected in the fact that the right side of ^ is un- 
bounded in any Sobolev topology. Hence we cannot get solutions of ([I]) by Picard 
iteration, and thus we arc not guaranteed smooth dependence on initial conditions. 

In Section|4]we study the dependence of a solution 77(1, s) on the initial velocity 
field w{s), given a fixed initial position 7(5). The Riemannian exponential map on 
is given by 

(9) exp^(w) = 77(1), where r] solves (H])-©. 

We will prove that for m ^ 3, as long as 7 e A™'^^, the exponential map is defined 
and continuous as a map from some open subset of T^A^ into A^. In fact we 
will show the exponential map is differentiable but not continuously differentiable. 
Our method is similar to that of Constantin-Kolev |CK| and Constantin-Kappeler- 
Kolev-Topalov jCKKTj : after establishing bounds on the linearized equation (to 
prove differentiability), we show that there are conjugate points arbitrarily close to 
by working out a very explicit special case (a string rotating like a rigid rod). 
If the exponential map were C^, then the fact that (dexp^)o is the identity would 
imply by the inverse function theorem that there is a neighborhood of on which 
there are no conjugate points. 

The failure of the exponential map to be has two consequences: one is that 
no geodesic can be minimizing, no matter how short (a conjugate point always 
implies the existence of a length-shortening variation); and the other is that we 
do not necessarily have geodesies joining two arcs (even nonminimizing geodesies), 
since the inverse function theorem is normally used to obtain this result. Hence 
the geometry of A"^ experiences some genuinely infinite-dimensional phenomena. 
On the other hand, the distance function generated is nondegenerate, since .4™ is 
a Riemannian submanifold of the geometrically flat space K."^. 

In Section [5] we compute the sectional curvature of .A™, showing that it is always 
positive. Intuitively, this implies stability of geodesies by the Rauch comparison 
theorem jCE| : however, the fact that it is unbounded above implies that the rigor- 
ous study of stability via curvature estimates faces some technical difficulties. In 
fact even if we could apply the Rauch theorem, the presence of conjugate points 
arbitrarily close to the identity makes it impossible to get any rigorous information 
about the growth of Jacobi fields. 

Finally in Section [S] we compare the geometry of A in the metric (jS]) to other ge- 
ometries on spaces of unparametrized curves, especially the metric on the space 
of parametrized curves modulo reparametrizations, studied by Michor-Mumford in 
jMMlj . Both metrics are too weak to preserve all the properties one expects in 
finite-dimensional geometry, but the metric ([8]) has a nondegenerate distance while 



(8) 




4 



STEPHEN C. PRESTON 



the Michor-Mumford metric gives a degenerate distance. We compare the geo- 
desic equation ([I])-© with the geodesic equation on A obtained from the Michor- 
Mumford metric: the primary difference is that our metric is essentially a sub- 
manifold metric, while the Michor-Mumford metric is essentially a Ricmmanian 
submersion metric on a homogeneous space. Since the submanifold metric ([8]) is 
related to the physical metric and has a nondegeneratc distance, we hope it may 
be of interest in shape recognition applications. 

In Appendix |X] we show how the results of this paper change if we consider 
periodic boundary conditions for the system (dl)-©. Many results actually become 
easier (for example, we can work in ordinary Sobolev spaces on the circle rather 
than weighted Sobolev spaces on the interval), and the essential features are the 
same. Then in Appendix |B] we explore what happens if we remove the constraint 
that all our odd curves have length 1; we see that many of the results break down 
in this case. 

Some of these results (in particular the nonnegativity of the sectional curvature 
in Theorem 15. ip were first obtained by Victor Yudovich, but not to my knowledge 
published. Alexander Shnirelman introduced me to this problem, and I thank him 
for many useful discussions on it. 

2. The manifold structure of the arc space A'^ 

For this section we assume all curves map into IR^, for simplicity. The space /C™ 
of odd curves in with the topology (|3]) is obviously a manifold, as a linear space. 
The topology defined by the seminorms (|4]) is sufficiently strong to make A™ a 
submanifold of /C™ when m ^ 4, but when m = 2 or m = 3 the topology is almost 
but not quite strong enough. The difficulty here is that a bound on the norms \\ri\\j 
for ^ J ^ 3 is not sufficient to ensure boundedness of sup_]^<j,<]^|77'(s)|, as shown 
by the example 

(10) r?'(s) = arctan (In (1 - s^)). 

First we recall the following lemma from |Prj . relating the weighted Sobolev 
norms ([5]) and the weighted supremum norm 

(11) lll/lll?fc= sup (1-.2).|/W(,)|2. 



Lemma 2.1. For any real j > and any nonnegative integer k, and any smooth 
function f, we have the following estimates for the norms ^ and 

(12) \\f\\U,k ^ ll/llL + 

(13) ml, < ii/iiL + 

Theorem 2.2. // m ^ 1, then the space defined by ^ is a C'^ Hilhert 

manifold. If m ^ 3, then A"^'^^ is a C'^' Hilhert submanifold of the space KT^^"^ 
defined by ([¥]). 

Proof If ry 6 /C^, then we have y_^{\ - s'^f\r^"{s)\'^ ds < CO. Hence by the standard 
Sobolev inequality on [— 1— e] for any e > 0, we see that rj is in [— 1-|-£, 1— e]. 
Hence it makes sense to impose the condition |r?'(s)| = 1 for s e (—1, 1), so that A^ 
is a closed subset of /C^. 

Furthermore for any such ry, we can write r]'{s) = ( cos 6* (s), sin 6'(s)), where 9 is 
uniquely determined once 9{0) is chosen. (Since 77' is continuous on (—1,1), so is 
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9.) We can easily compute that we have r] e if and only if e J-'^[M], where 

(14) J-^P] = |/: [-1,1] : J (1 - s2)J+i|/(j)(s)p ds < oo for ^ j ^ mj 

for F = M or C. Now J^^ is a Hilbert space, and the map rj ^ 9 defines coordinate 
charts (for example, on the set of rj with 7/(0) ^ v for any fixed unit vector w), 
for which the coordinate transition maps are trivially C"^'. In this way we get a 
smooth Hilbert manifold structure on A? , and the same process will give a manifold 
structure on any A™"^^: we just have to check that r] 6 ^™+^ if and only \i 9 e 
which will follow from the next result. 
To actually obtain ^™+^ as a submanifold of /C^^^, we can construct a co- 
ordinate chart on /C™''"^ which makes this obvious. We simply write r]'{s) = 
(e'^(^) cos6(s),e'^'^^^ sm9{s)) for functions -0 and 6. Obviously if we think of 
as C, this is just the exponential ri'(s) = e^'"' where ^(s) = ■i/;(s) + i9{s). Now our 
claim is that if ip is bounded (equivalently, if \ri'\ is bounded both above and below 
away from zero), then rj e /C™ if and only if ^ e ^™~^[C]. This is relatively easy to 
check using the Faa di Bruno formula for the derivative of a composition: we have 
T]' = and ^ = In 77' , so that for m ^ 1 the formula yields 

(15) ^-e^y^^nr^^^'^'' 



7- 1 



where in the second sum we set £ = ki + ■ ■ ■ + fc,„ , and where both sums are taken 
over all nonnegative integers kj such that I ■ ki + ■ ■ ■ + m ■ km = m. 

We can now prove for ttt, ^ 1 that ^ e J^'"[C] and supSR^ < 00 implies 77 e /C'"^^. 
Ignoring the specific constants and using Cauchy-Schwarz, (|15p yields 



l'7llm + l,m + l ~ 



pi m 

supe2^«;^ (1 - []|C(^")(s)l''^ ds 

J-i 

supe^^«2 I (l--')n((l-^')^l^^^'(^)l') ^ ds. 



Now if j ^771 — 2, we can pull out the term |(1 — s^)^|^^^'p, using Lemma [2.11 to 
get 

III^IIIL- ^ II^IIL + + l ^ 11^11?+!,. + ll^ll?+2,, + l + IICII^3,,+2 < 00 

since j + 2 m and ^ e J-"™ . Hence we only need to worry about bounding 



(17) J (1 -s2)" + l|^("-l)(s)|2'="-i|^(™)(s)|2'c,„ 

Since {m — l)km-i + inkm ^ ttt,, we cannot have both km-i and km nonzero. It 
is then easy to check that in all possible cases, the term (|17p can be bounded in 
terms of the J^'" norm of ^. 

We can similarly prove for 771 ^ 1 that rj e /C^+i and inf jTy'l > implies ^ e J"™, 
by using formula (|16p . We note as a consequence that when ^ is purely imaginary 
(corresponding to 7; e A™), we get a bijective correspondence between an open 
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subset of y^™"*"^ and an open subset of {0} x for any m ^ 1. This yields 

coordinate charts on ^™+^ for m ^ 1. 

Now to show that ^"'+^ is a submanifold of /C^+i for to ^ 3, we want to show 
that TO ^ 3 imphes that the coordinate chart rj' gives a bijection between 

^m+i g^j^j J^™[C], and the only thing remaining after the above estimates is getting 
upper and lower bounds on |?/| and 5R^. Ideally we would use p^ . but that estimate 
fails since our case corresponds to j = (as mentioned in (|10p ). Instead we use the 
usual (unweighted) Sobolev inequality on [—1, 1] to get 



sup i^(s)p < C ms)?+ C w{s)\^ds 

-l=Sss=l J-1 J-1 



\l.o + \mii, 



and then apply ([T2l) to get 



IllV'llllo^ll^lllo + IIV'llli + IIV'llf.i + llV' 

(18) 



< 



V^IIL + Il^llli + II^IIL + Il^f 



11,0 ^ II'/^II2,1 ^ II'/^II3,2 ^ Ilrll4,3- 

Thus if i/j = 5R^, the fact that ^ e J^'^[C] implies ip is bounded. Similarly \ri'\ can be 
bounded in terms of the norm of /C*, so that the space of curves with \ri'\ bounded 
away from zero is an open subset of /C^. 

We conclude that for to ^ 3, the map r/ ^ is a coordinate chart on /C™"*"^ which 
models /C™+^ on the Hilbert space J^'"[C]. Furthermore locally the set ^™+^ is 
mapped bijectively under this coordinate chart to the closed subspace {0} x J^"'[R], 
so that ^"'+1 is a Hilbert submanifold of /C™+^ □ 

We note that on the way, we actually establish that ^™+^ is a Banach subman- 
ifold of the slightly modified space IC"^^ for to ^ 1, where the topology on /C™'*' 
is given by the Banach norm 



SIM 



j=a 

For TO ^ 3 the spaces /C™"*"^ and K™^^ are isomorphic since the extra term |||77|||o i 
becomes redundant due to (|18p . 

We can show directly that A"^ is not a Hilbert submanifold of /C^ (and similarly 
that A'^ is not a Hilbert submanifold of K,^) by the following calculation: we give A^ 
the coordinate chart J-^ and view /C^ as a Hilbert space with the identity coordinate 
chart, and then the embedding of A^ into /C^ is given in coordinates by the map 
2: 9 1-^ (cos 0, sin 0). The derivative of I at w e TqJ^^ = J^^ must be given by 

C = DX{e) (cj) = ( - cj sin 61, w cos 6) . 

However if a; e J^^ and 6* e J^^ , we do not necessarily have C 6 /C^ , since 

IICII^,2 = {l-Sy{u:{sf+J{sf+Uj{sfe'{sf) ds, 

and we need a bound on sup a; in order for the last term to be bounded by ||0||2 i- 
But there is no reason any lo € has to have bounded supremum, due to examples 
like (fTOj) . Hence the embedding of A? into K? is not even differentiable, and so A? 
cannot be a smooth submanifold of . 
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3. The geodesic equation 

Since the space /C™ defined by (jlj is a Hilbcrt space, it has an obvious Rie- 
mannian metric given by the Hilbert norm. Geodesies in this metric are always of 
the form ri{t) = 77(0) + r]'(0)t, so the exponential map is exp^(w) = 7 + tw. This 
formula is the same regardless of whether we define the Riemannian metric by the 
weighted Sobolev norm ([5]), as in 

~1 m 

(19) «u,«>X_„= Y,(^-sy(D^u{s),D^v{s)yds, 

or the weaker norm given by the kinetic energy, 

(20) {(u,v})^ = j (u{s),v{s)}ds, 

which of course is (|19p when m = 0. However, geodesies on the submanifold A™ 
will be different depending on which choice we make. The natural choice from the 
perspective of physics is ^U\\ . which is a weak metric (it is not equivalent to the 
Hilbert norm). 

Formula ([20| obviously gives a smooth metric on the manifold JC"\ while as 
discussed in Section [2] the space .A™ is a smooth submanifold of AC™ if m ^ 4; 
nonetheless the connection on is not smooth. The reason for this is that the 
connection of a submanifold is obtained from the connection on the full manifold by 
orthogonal projection, and the orthogonal projection operator is not smooth (using 
the metric ((20|)). In other situations when we get a smooth ODE on a submanifold, 
it is due to smoothness of the orthogonal projection: see for example 'EM' and 
pP] . 

3.1. The orthogonal projection. We will show in this subsection that the or- 
thogonal projection is intimately related to the ODE ([2]). Since it will appear re- 
peatedly, we summarize the main properties of the Green function for it, as proved 
in [P?] . 

Theorem 3.1. Let G: [—1,1] x [—1,1] K denote the Green function of 
satisfying 

(21) ^-\y'(s)\^Gis,x) = -5{s-x), G(l,x)=G(-l,x)=0, 
so that the solution of ([2]) is 

(22) ^(s) = J G{s,x)\Tjt^fdx. 

Then G{s, x) = G{x, s) for all s and x. In addition G(s, x) ^ for all s and x, and 
G{s,x) = only on the boundary. Furthermore if G{s,x) = ^^G{s,x) + G(s, — x)] 
is the even symmetrized Green function, then 

(23) G(s, x) ^ min{l - |s|,l - 
and 

(24) G{s,x)^{l-\s\){l-\x\)e-y{l + g) where g = { {1 - IsDltjssf ds. 

Jo 
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Proposition 3.2. Suppose m ^ 2. The orthogonal projection Pry : T^K™ T^A^ 
generated by the metric (|20p is given by 

(25) P-yi^) = ^ ~ i'^l'y where a solves 

(26) a"{s) - |7"(s)| V(s) = <z'(s), 7'(s)> with a(-l) = a(l) = 0. 

For any fixed 7 6 A™'^'^, the projection is continuous from T~^IC™ to T-yA™ , but 
the map is not continuous from TIO^ to TA"^ . 



Proof. Wc first observe that for any function a with (t(1) = o-(— 1) = 0, the vector 



field V = -r-(o'7') is orthogonal to T^A"^ in the metric (PO)) . since if (w',"f') = 



then 

((v,w)\ = ^(w{s),f^{<j{s)^'{s))}ds 

= {a{s)(w{s),^'{s)))^='_^ - ^a{sKw'{s),j'{s))ds 
= 0. 

Now given z e T^/C™, solve the ordinary differential equation ([26]) for a, and 
define w = z — {a^')' . Then we have 

{w'{s),i{s)) = <2'(s), y(s)> - <a"(s)7'(s) + 2a'{s)Y{s) + ais)Y" {s),j' (s)) 

= (z'is),y{s))-a"{s)+a{s)\Y{s)f=0, 

using the fact that |7'(s)P = 1 implies that (^'{s), 7"(s)) = and hence (7'(s), 7"'(s)) 
— |7"(s)p. So w actually satisfies the tangent condition (w',Yy = 0. 

We just need to check that w actually is in TjA^, i.e., that < 00 for 

1 ^ j ^ m as long as 2; e T^KT^ for m ^ 4. It is obviously sufficient (and easier) 
to check that v e T^KT^ . Our estimates are based on the same estimates that are 
done in |Prj . to which we refer for more details. 

The key is that, by the product rule, v = {cr^')' satisfies 

m + l 



m + 1 pi 
I — n -^ — 1 



(27) 



^ ^^(ll7llm + 2,m+2 + Il7llm,m + l) + Ill7 III m-^-l,m-« Ik 



2 

£=0 



where A = sup^lcr'l ^ supj,|tT|/(l — s^). Based on the bound for the Green 
function for (pS)). we easily see that 

^'^lkllo,l^Pllll + Pll2,2, 

using ([T2|) . We easily get 

lll7lllm-£-l,m-f ~ ll7llm-£,m-« + ll7llm-£+l,m-f+l + ll7llm-£+2,m-<'+2 

using ()12p -(fT3 |) . which allows us to bound ||w||m „ in terms of the ^'"+'^ norm of 
7, once we get a bound on ||cr||^_^^ ^_^2 for ^ ^ ^ m — 1. 

To obtain this, we use ((26|) along with the product rule and Lemma [2.11 to get 
the recursive inequality 

ikii'+M+2 ^ (sfi'ii7iiy'(2:-.2ikii?-i,, +s-;;ikiiy. 
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Combining this with ([27|) gives 




m /m+2 
i=0 \i=l 

for m ^ 2, for some function Gm- 

On the other hand, if 7 6 and not ^™+^, we obviously do not in general 
have {crj'y 6 T^/C™ even if a is C^ . Hence even if z is C^ , the projection -P-y(2;) 
given by ((25)) is not in T^A™' if 7 is only in A"^ . □ 

3.2. The second fundamental form. The orthogonal projection encodes all the 
geometry of the submanifold A™" , via the second fundamental form. The second 
fundamental form then leads to both the geodesic equation and to the sectional 
curvature (which we will discuss in Section [S]) . 



Lemma 3.3. The second fundamental form of A"" , as a submanifold 0/ /C™, is 



given by the operator S: T^A"'' x T^A"" (T^A^y defined by 



(28) S{u,v) = -^{(Juvis)j'is)), where 



<Juv{s) = J 



(29) auv{s)=\ G{s,x){u'{x),v'{x)}dx 



in terms of the Green function of Theorem \3.1\ This operator is only well-defined 

Proof. In general the second fundamental form is the orthogonal projection of the 
connection: if U and V are vector fields on A"^, with u = and v = the 
values of these fields at 7 e A""", then the second fundamental form is S{u,v) = 
(^(\/ijV)-yj , and the value obtained depends only on the values U-y and (not on 
the extensions U and V). (See for example do Carmo |dC2j .) Unfortunately it is 
somewhat awkward to work with general vector fields on a function space, and so 
we use the alternative method of vector fields along curves. 

So suppose ri{t) is a curve in A"^ with 77(0) = 7, and let V{t) be a curve along ?y, 
so that V{t) 6 Tj^(t)^™ for each t. Set u = 77(0) and v = V{0). Then the covariant 
derivative of V in the direction u, calculated in the fiat ambient manifold /C™, is 

where the last equality comes from using fiatness of the Hilbert manifold /C™ to 
change the covariant derivative to an ordinary derivative. Now using formula (j25p 
for the orthogonal projection, we get 

cT'Us)-\j"{s)\'=(Vts,Vsyt=o- 

To simplify (T^si'ys), we use the fact that V{t) e T^(t)yl™ at each time to obtain 
(Fs, ?7s) = for all time, so that diff'erentiating we get (Vst, ils) + (^s, 'Hst) = 0. Now 
at time we have Vs = v' and rjst = u' , so that (Vts, 7/^)4=0 = —(u',v')>, which 
implies 

oiM-\i\^)?^u.{s) = -{u\s),v\s)y. 

Hence the formula ([22]) yields □ 
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As a consequence we obtain the geodesic equation on A™, using the general 
formula for a geodesic on a submanifold: 

which using ((28)) reduces to ^ with tension given by ([2]). 

The fact that the orthogonal projection fails to be continuous in both 7 and 
z implies that, unlike in Ebin-Marsden |EM| . the weak geodesic equation on A"^ 
cannot be viewed as an ordinary differential equation, and hence cannot be solved 
via Picard iteration. 

On the other hand, we can prove local existence and uniqueness of solutions. The 
main result of the author's companion paper |Prj is the following theorem (restated 
here more geometrically): 

Theorem 3.4. Suppose ni ^ 3. If ^ e ^™+^ and w e T^A'"', then there is a 
T > such that there is a unique solution r): (— T, T) —>■ ^'"+^ of ([I])-© satisfying 
rj{0) = 7, 77(0) = w, and such that r]{t) 6 A"'+'^ and T)(t) 6 T^(t)A'''- for all t. 

The loss of derivatives here (i.e., the fact that 77 is not as smooth as t]) means 
that we do not have a classical exponential map (which would be a map from an 
open subset of TA"^ to itself) . 

However, if we fix an initial configuration 7 e A™''^^ , then wc have a reduced 
exponential map 

(30) exp^ : Vl c T^A" A"" 

defined on some open neighborhood i7 of by exp^{w) = 77(1), where 77 solves ([1])- 
([2]) with 77(0) = 7 and 77(0) = w. Actually as mentioned exp^(w;) is really in A"^'^^, 
but the theorem in the next section on continuous dependence can only be proved 
in this weaker topology. 

4. Differentiability of the reduced exponential map 

We now want to establish continuity of the reduced exponential map pO|) : in 
other words, for a fixed, sufficiently smooth initial whip configuration 7, we show 
that the solution depends continuously on the initial velocity w in any weighted 
Sobolev topology ([5]). In fact we will obtain Lipschitz continuity as a result of 
showing that exp^ is differentiable everywhere on T^A"^, with bounded derivative, 
but that it is not C^. The latter fact is a consequence of clustering of conjugate 
points near 0, the same thing that happens for the exponential maps corresponding 
to Burgers' equation jCKj and the Korteweg-dcVrics equation jCKKT] . 

First we compute the derivative of the exponential map (which is just the lin- 
earization of the equations ([T])-©) and prove that it is bounded. 

Theorem 4.1. Suppose m ^ S and 7 6 ^"'+^, and let w 6 T^A"^ . 

The derivative 0/ exp^ : T^A"^ —>■ A"^ is D{cxp^)w{y) = ^(1) where i^,ip) is a 
solution of the linearized equations 

(31) (tt = dsiai.) + dsiM. 

(32) (l)ss - |77ssP'/' = 2{iiss,^ss}o- - 2(rist,^st) 

where (^CsiVs") = 0, with initial conditions ^(0, s) = and ^t(0, s) = y{s). Here rj 
and a solve with initial conditions 77(0, s) = 7(5) and 7/t(0, s) = w{s). 
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The derivative satisfies the bound 

ra 

(33) 2 P(exp,).(y)|ll. < Urn (SSl7lL?„E,"lilkllL-) Sllolbll^ 

fc=0 

/or some junction Hm- 

Proof. Wc obtain (|5T|) - ([5^ by considering a family of solutions (?7(r, t, s),a{r, t, s)) 
depending on a parameter r, satisfying 77(r, 0, s) = 7(5), 774(0, 0,s) = i«(s), and 
r]tr{0,0,s) = y{s). Setting ^ = 77r|^_Q and (j) = <^r\j._Q, we get the desired equations 
and initial conditions. 

The bound (j33p will be obtained by bounding the tension- weighted energy norms 

m — 1 ^1 

(34) em-i=2 (a{t,sY\didtat,s)\^ +a{t,sY+'\di+'at,s)f)ds. 

As in the estimates of jPrj . we compare the tension- weighted norm to the time- 
independent weighted energy norm 

m — 1 ^1 

e=o -^-1 

The bounds from Theorem 13 . II imply, as in jPrj . that we have bounds A{t) and B{t) 
such that 

1 cr(t, s) , , , 
0<^^A^^A(t)<co 

implies that the norms ([M)) and ([55)1 are equivalent. We recall also the result from 
[Pr] that supg|(Tt(i, s)|/(l — s^) ^ supg|crst| = C'(i) is bounded in terms of the 
norm of 7 and the norm of w. Finally we write 
m ^1 

Emit) = Y, ({l-sy\d%\' + il-sy+'\8i+'r^\Ads. 

By our assumption on the initial conditions and Theorem 13.41 Em{t) is bounded. 

Now we compute the time derivative dim-i/dt one term at a time, obtaining 
after some simplifications 



dt _ 

i+i 



s^Y\di<j\^\dl+^-'^\^ds 



J=2 
i+1 



+ A'Uthi 2 /I - s2)^|5^0P|5f+'-^r/P ds 

+ A'^A\\^t\U,e{\\vh+2,e+2 + hll«-n) 

+ f^ds{a'^\dl+'C,d%))ds, 

where A = sups|0s|. 

The terms on the first line are obviously bounded by £m-i, and the terms on 
the last line vanish because of the boundary condition on a. The fourth line is 
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bounded in terms of Em and Em-i, using a bound on A obtained by using the 
Green function bound ([23|) . together with the weighted Sobolev bound IE]), to get 



< J ^(1 - s^)\S,ss?\'nss? rfs + J \U\'\Vst\'ds < {A' + l)e2E2. 

The terms on the second hue can be bounded using the fact from [Pr] that ||cr||^_,_i 
can be bounded in terms of Ej for j ^ 3. FinaUy the terms on the third hue can 
be bounded by obtaining a bound on Hi/ill^.^ in terms of e^-i, which is obtained 
in the same way as the proof of Proposition 13.21 

Summing from £ = to to — 1, we obtain an estimate of the form 

— $ J{rjmj£m-l 

for some function J, and then the bound 1 — ^ Ba gives us a Gronwah inequahty 
of the form dira-i/dt ^ J{Em)£ra-i- We then obtain the desired bound ([33]) from 
this. □ 

Integrating the derivative obviously gives us a bound on ||exp^(w)— cxp^(i;)||m 
as in Cheegcr-Ebin jCEj . which estabhshes that the reduced exponential map ([30]) 
is locally Lipschitz, as desired. 

However, the exponential map cannot be C^; if it were, then the fact that its 
differential is invertible at zero would imply it is also invertible near zero. But the 
latter does not happen. 

Theorem 4.2. The reduced exponential map (j30p is not a map on T-yA"^ for 
any m. 

Proof. For any 7, the differential (Dexp^)iL, at = is the identity, which is easy 
to see from the fact that in this case, ri{t) = 7 and (j{t) = for all t in ([T])-®. This 
implies that the solution of ([32]) is zero, so that ([3T|) reduces to £,tt = 0. Since 
^(0) = and ^(0) = y, we get ^(1) = y, i.e., (Dexp^)o(y) = y. If the exponential 
map were C^, then (Dexp^)^, would have to also be invertible for sufficiently small 
w. 

However, we can find w arbitrarily close to in T^A"^ for which (ZJexp^)^, is 
not an isomorphism. To do this, we work out an explicit solution in detail. It is 
easy to verify that 

„2\ 



s 



(36) r]{t,s) = {s cos ojt,s sin ojt) and a-{t,s) = ^{1 
form a solution of ([T])-® for any angular velocity w. In this case we have of course 

(37) 7(s) = (s,0) and w;(s) = (0,ws). 

The constraint <(^s,77s) = implies that ^^(t, s) = ■s)(~ sinojt, coswt) for 
some function x- Differentiating pip with respect to s gives (j)ss = and 

(38) ^^^.^\ = ^3U{l-s')x). 

Since ^ is odd as a function of s, wc must have x even as a function of s. Expanding 
X in a basis of derivatives of odd Lcgendre polynomials as 



X{t,s) = £ Xn{t)PL-l{s) 
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■ ((1 - s^)PL-i{s)) = -2n{2n - l)PL-i{sl 



and using the fact that 
we see that 

X':,{t) - u:\n{t) = -UJ^n{2n - l)xn{t), 
the solution of which, with Xn(0) = and Xn(0) = c„, is 

Xn{t) = — sin(a„t), 

an 

where = uj^{2n + l)(n - 1). 

So if y„(s) = (0, F^„_i(s)) for any n, then we have (D exp^)^(?;„) = ^^^yn- 
Now for each ?i ^ 2 we can define a;„ = -y/7r/((2n + l)(n — 1)), so that = tt. 
Obviously a;„ ^ and so the corresponding w„ in ([57]) converge to in TjA™' 
for any m, yet (Dexp^)^^ has a nontrivial kernel for any n ^ 2. Hence even for 
a fixed smooth 7, the differential w i—>- (Dcxp^)„ cannot be continuous as a map 
from T^A"" to L(T^^", T^™). □ 

Generally speaking, if t] : [0, T] M is a geodesic in a finite-dimensional mani- 
fold M such that 77(6) is conjugate to 77(0) for some b ^T, then 77 cannot be locally 
minimizing on [0,r]. (See for example do Carmo [dC2] .l For the geodesic ([M]). 
our computation shows that no matter how small T is, there is a conjugate point 
at some 6 < T, so that even an arbitrarily short geodesic cannot be minimizing. 

On the other hand, the induced distance is not degenerate: if 71 and 72 are 
distinct curves in A, then the infimum of lengths of paths in A joining them has 
a positive lower bound. This shows that the vanishing geodesic distance is not a 
consequence of having unbounded curvature, which was sugg ested in |MM3| . 

Proposition 4.3. The Riemannian distance between distinct curves 71 and 72 in 
A, in the metric (jSj, is always positive. 

Proof. The idea is basically that the Riemannian distance in a submanifold of a flat 
space is always at least as large as the "chord" distance in the flat space. Specifically, 
if r](t,s) is a curve with \r]s\ = 1 and ry(0, s) = 71 (s) and 77(1,5) = 72(5), then the 
length of 77 is bounded using the Cauchy-Schwarz inequality by 

r]t(t, s) \ dsdt 
\l2{s) -ji{s)\ds. 

So we get an absolute positive lower bound regardless of 77. □ 
5. Curvature of the arc space 

Having computed the second fundamental form of A™ in /C" in Lemma 13.31 we 
can now compute the sectional curvature using the Gauss-Codazzi formula. 

Theorem 5.1. Let m ^ 2 and 7 6 ^"'+2, and let u,v e T^A™ be vector fields. 
Then the sectional curvature K{u, v) at 7 in the section spanned by u and v is 

l'_^l^_^G{s,x)(\u'(s)\^\v'{x)\^ -{u'{x),v'{x)Xu'{s),v'(s))dsdx 




(39) K 



1^11*^(^)1^ 1^(2^) P ~ '\'U'ix),v{x)')(u{s), v(s)yds dx 
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where G is the Green function (|2ip . The curvature satisfies K ^ e ^/(l + g), where 
g = ^^(1 — s)|7"(s)p ds; in particular it is always positive but never bounded above. 

Proof. First wc note that /C™ is flat, so the Gauss-Codazzi formula gives 
((S{u, u), S{v, v)}} - ((S{u, v), S{u, v)}} 
((u,u}X(v,v)}-((u,vy 
where S is the second fundamental form (1281). Now we have 



((S{u, u), S{v, v))) = fWuu^ + ^uul\ a;„7' + ^v.i')ds 

= J ^'uu^'vv + \l"\^(^uu(Jvv = J <yuu{s)\v' {s)\^ ds. 



Now by (HSl) we have 

o'uu(s) = J G{s,x)\u'{x)\^ dx, 
where G is the Green function (PT|) . Hence we can write 

P ds dx. 



{{S{u,u\S{v,v))) = ^^ ^ G{s,x)\u'{x)\^\v'{s)\^ 

Formula ([55)1 follows. 

Nonnegativity of the sectional curvature follows from the fact that, by symmetry 
of the Green function, we can write the numerator of (1391) as 



-J" J" G{s,x)M{s,x) ds dx, where 



2 . _ . 

M(s, a;) = \u'is)f\v'{x)f + \u' (x)f \v' is)f - 2(u'(s), w'(s)><u'(x), v'(x)). 

We have G{s,x) ^ for all s and x by Theorem 13. II and M{s,x) ^ for all s and 
X by the Cauchy-Schwarz inequality. 

To get a sharper estimate, note that since 7 e by assumption, we know ||7||i,2 
is bounded by ([T2|) . Hence we can use the estimate ([24]) to obtain 



on [—1,1] X [—1,1], which allows us to write (|5^ as 

^ ' - Ail + g) \\u\\lo\\v\\l,-((u,v})l. 

To get a lower bound on this, we disregard the restrictions on u and v (that they 
be elements of T^A'") and minimize over all odd vector fields along 7. Expand 
u and V in a basis of odd Legendre polynomials as u{s) = Xin^i '"n^2n-i(s) and 
^'(*) = 2n=i ^n^2n-i(s), whcre u„ and Vn are vectors in M.^ . Then the bound ([40)1 
becomes 

e"'? E^.m^i A„A,„(|u„p|u„ip - <u„,-y„XMm,'ym>) Afe"'? 

^ — ^ ; '-^ T. — — ; ^ 



4(1 + q) i]::„=i (I 

p - <w„,w„><u„,'!;„» 4(1 + 

where A„ = 2n(2n — 1). Positivity of the curvature follows. It is easy to see that 
the curvature can be made arbitrarily large using this formula as well. □ 
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The fact that the curvature is unbounded above is responsible for the fact (as 
shown in Theorem 14. 2 p that conjugate points afong a geodesic occur at times arbi- 
trarily close to 0, and hence for the failure of the Riemannian exponential map to 
be C^: If the curvature were bounded above, then the Ranch comparison theorem 
would imply that there is a small interval of any geodesic in which no conjugate 
points can occur, contradicting Theorem 14.21 



The space of curves is of interest in shape- recognition applications |MMlj , since 
the first step in distinguishing two shapes in the plane is to distinguish their bound- 
ary curves. Obviously in studying geometry of curves for this purpose, we want to 
consider the image of the curve in the plane (which is all the viewer can see), not 
the actual map from the interval to the plane. There are essentially two ways to 
get a Riemannian structure on this set: impose a specific parametrization with unit 
speed (as we have done in this article so far), or consider all parametrized curves 
and mod out by the reparametrizations (the diffeomorphism group of the interval). 
The latter approach is the one taken by Michor and Mumford |MM1| . However the 
approaches are basically equivalent an odd curve of length 2 in the plane always 
has exactly one parametrization on [—1, 1] of unit speed, so that if 



Imm2([-1,1],M2) = {r/: [-1,1] ^R^ : |r/(s)| #0, ry(-s) = -ij{s), i(r/) = 2} 



denotes the space of odd immersions into for which the image has length 2 and 
2?([— 1, 1]) is the group of odd orientation-preserving diffeomorphisms of [—1, 1] to 
itself, then we expect to have Imm2([— 1, 1])/X'([— 1, 1]) = A. This doesn't quite 
work rigorously since the action is not always free, but we can still see what the 
Michor-Mumford metric looks like on A. In this section we will assume all objects 
are C'^' and work formally, although with a bit more work we could extend the 
results to weighted Sobolev spaces. 

The reparametrization action of !?([— 1, 1]) on Imm2([— 1, 1]) is given by compo- 
sition: for h 6 !?([— 1, 1]) the map is Rhiv) = o h. Wc can define a Riemannian 
metric on Imm2 by 



We clearly have ((u oh,vo h^rioh = {(u, v^tj for any e !?([— 1, 1]) by the change of 
variables formula, so that the metric (j42p is invariant under the group action. This 
is in contrast to the metric ([5]), which is not invariant under reparametrizations. 
(Of course on the submanifold of unit-speed curves A, both metrics coincide.) The 
geodesic equation in the metric ([5]) is just r]tt = 0, while the equation in the invariant 



Here we modify the Michor-Mumford space to consider curves with a free boundary and a 
fixed length; in Appendix|A]we consider the periodic case, and in Appendix[B]we discuss the case 
where length is not constrained. 



6. Comparison with other metrics 



(41) 



(42) 
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metric ()42p is the much more comphcatecd nonhiiear eUiptic equation 

6.1. The Michor-Mumford metric on the arc space. Now we formally iden- 
tify the quotient Imm2 I'D with the space A in order to compare the induced metric 
on A to our metric ^ . 

Theorem 6.1. Define the standard reparametrization map $: Imm2 A by 
$(77) = 11 o where h 6 2?([— 1, 1]) is given by h{s) = ^^Irj' {x)\ dx . Then for 
any k 6 1, 1]), we have $(77 o k) = ^{rf), so that $ is invariant under the 

action ofD. Hence descends to a map from the quotient space ImmijD into A. 

There is a unique metric on A defined by the condition that $ is a Riemannian 
submersion, and it is given by 

(43) = ^\u{s),-i\s)^){v{s)n'{s)^)ds, 

where u and v are in T^A and ^'{s)-^ is the rotation of the unit vector Y{s) in 
by 90°. 

Proof First, if 7 = $(77) then y{s) = , so that |7'(s)| = 1. Hence $ 

actually maps into A. Now if k: [—1, 1] [—1, 1] is a difFeomorphism with fc' > 0, 
then for the curve 770A: wc get h = hok,so that 7 = (770^)0 (hok)~^ = rioh~^ = 7. 
Hence (f>(7y o k) = $(77), so that $ is invariant under the reparametrization action. 

To get the metric, wc first compute the derivative Z3<f>. For any odd vector field w 
along an odd immersion 77, let s) = 77(5) +ew{s); then for sufficiently small e the 
map s i—> x(e, s) is still an odd immersion, and we have (D$)^(7«) = ^|^_p<f>(x(^))- 
It is then easy to compute that 

(44) {D^\{w) = /3 o /i-i, where 

Pis) = w{s) - {i;(w'{x),r,'{x)/Wix)\}dx)v'{s)/Wis)\. 

If 7 = $(7/) = $ o h~^, then wc can check that z = (13$),, (tjj) actually satisfies 
(z', 7') = as expected for any w, and that the kernel of (£'$)r; is the vertical space 

V„ = {/77':/: [-1,1]^ Mis odd}. 

The horizontal space is the orthogonal complement of the vertical space in the 
metric (|42|). which is 

H„ = {f{vY-f- [-1,1]^ Mis odd}. 
The metric on A which makes $ a submersion is given for z 6 T-yA by ((z, z)).^ = 
<^w, w'^ri, where 77 is any curve with $(77) = 7, w e T,,Imm2 is any horizontal vector 
field with (D^)ri{w) = z, and the right side is computed using the invariant metric 

^This is a typical sort of tradeoff for invariance: the same tiling happens in fluid mechanics, 
when we consider the diffeomorphism group ©(M) of a Riemannian manifold M and the volu- 
morphism group ©^(A/) = {r} e 'D(M) : rf* = /i} where ^ is the Riemannian volume form on 
M. The simplest metric on "D^M) is the non-invariant metric {(it o 77, i; o rf^r^ — ^j^^UyV") o rjdfi, 
for which the geodesic equation is ritt = (which leads to Burgers' equation ut + V^f = 0, where 
rit — u o rf). The right-invariant metric is {(it orj,v o rjj^r) ~ d-fj,, on which the geodesic 

equation is itt -I- \'uU + (div «)u -I- iVjup — 0, again with rit — u o rj. Both metrics agree on the 
submanifold V^IM), on which the projected geodesic equation is ut +\'uU = — Vp with div it = 0. 
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(|42| . Invariance of the metric (|42|) ensures that we get the same {(z, z))-^ no matter 
which 77 we use, so we might as well use 77 = 7. Then wc can compute that the 
unique horizontal w with (D^)^{'w) = z is w{s) = (z(s), 7'(s)-'-)7'(s)-'-, and formula 
(I4l follows. □ 



Next let us compute the geodesic equation on A in the Michor-Mumford metric. 
The following lemma is helpful in finding compatibility conditions for it. 

Lemma 6.2. Ifij: [0,T] x [—1, 1] — > is a smooth curve, and if we define £ = \r]s\ 
and 

then we have the compatibility equations 

(46) if = bs — aK£, dt{£K) = ujs, md Og = £u) — bK£ 

Proof. Wc write 



Differentiating rjt with respect to s and matching coefficients with rj^t, we get as/£ = 
uj — bn and £t = bg — aK£. Then using rjsts = Vsst we obtain dt{K£) = t^s- D 

Theorem 6.3. A geodesic rj in A with the metric (j43p satisfies the equations 

(47) at = ^Ko^ + btts, Kt = ojs, bs = na, Ug = oj — nb, 
where 

(48) a = (iit,V^}, b = (r]t,Vs}, k = (Vss,'ns}^ ^ = (Vst, Vs)^ 
and a(— 1) = a(l) = 0. 

Proof. If rj is a curve in A and u is a variation field along rj, then it is easy to 
compute that the first variation of energy in the direction u is 

(49) «'"^i«'7t»^=!_i + {u^,atris - asrit}dsdt. 
Jo Jo J-i 

For 7y to be a geodesic, this must vanish for every u e T^,4, i.e., whenever (u, 77') = 0. 
From the first term we get the boundary condition a(— 1) = a(l) = 0, and from 
the second term we get the equation at = bag + ^na^. The other equations are 
obtained by setting £ = 1 in Lemma [6.21 □ 

The geodesic equations take a slightly different form than that given in |MMlj : 
there the authors use the normalization 6 = rather than our normalization £ = 1. 
Of course, the images of the curves in are necessarily the same. 

The drawback of the Michor-Mumford metric, as discussed in |MM1| . is that 
the induced Ricmannian distance between elements of A is zero; that is, for any 
pair of curves 71 and 72 in A and any e > 0, there is a curve rj in A with 77(0) = 71 
and 7;(1) = 72 such that II '^^ ^ ^- shown in Proposition 14.31 our metric 

on A does give a genuine nondegenerate distance. 

We now ask what a right-invariant metric on Imm2([— 1, 1],IR^) would have to 
look like in order to give our metric ^ on A as a, Riemannian submersion using 
the procedure in Theorem 16. II 



£. 
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Theorem 6.4. Let us define a Riemannian metric on /mm2([— 1, 1], K^) as follows: 
if 1] is a curve and w is a vector field along rj, let 

(50) ^^(.y '^-+Li'^^(-)iUo iv(x)i 'V 

Then the metric (jSOp is invariant under the reparametrization action by 2?([— 1, 1]), 
and the map Imm2 A defined as in Theorem \6.1\ is a Riemannian submersion 
onto the arc space A in the metric ([8]). 

Proof. To check invariance, we just need to verify (^w ok,w o fc)),;o/j = ^w, w)),, for 
any increasing diffcomorphism fc of [— 1, 1]. This is straightforward from the change 
of variables formula. 

To check the submersion condition, we suppose we have a curve 7 with |7'(s)| = 1 
and that w is a horizontal vector field along 7, i.e., that <(w(s), 7'(s)) = 0. Then as 
in Theorem 16.11 we have 

{D<^Uw){s) = w{s) - {%{w'{xli{x))dx)i{s). 

Suppose z 6 T^A, i.e., that z is a vector field along 7 with (z'(s), 7'(s)) = 0. Then 
z{s) = f{s)'-f'{s) + g{s)y{s)-^ where f'{s) = K{s)g{s), and to get {D^){w) = z 
where w is horizontal, wc must have w{s) = g{s)'y'{s)-^. We can then check that 
the definition ([5(1)) yields 

((w,w}\= j'^[f{sf+g{sf] ds, 

as desired. □ 

Of course there are other choices for (|50p ; only the inner product of horizontal 
vectors is determined by the submersion condition, and wc can use any formula at 
all for vertical vectors. 

6.2. The metric on A. Finally we relate both our metric ([8]) and the Michor- 
Mumford metric ([43]) to another choice of distance on the arc space A. Klassen et 
al. |KSMJj pointed out that unit-speed curves in can most easily be represented 
in terms of their angular representation 6 defined by 7'(s) = ( cos 6'(s), sin 0(s)), as 
we did in Theorem 12.21 to get a coordinate chart on A™. Since 7(0) = in our 
space, we obtain 7 by integrating: 



(51) 7(s) 



cos9{x) dx, j sui0{x)dx 



Since the space of (even) functions 9: [—1, 1] ^ M is a linear space, it has a simple 
choice of Riemannian metric arising from the standard Hilbcrt structure. That is, 
if 7 is a curve with angular representation 0, and w is a vector field along 6, then 
the KSMJ metric is 



(52) «w,w» = J uj{sfds. 

In the physical space E^, the angular tangent vector w corresponds to the vector 
field 



u(s) = ijj{x)Y{x)'^ 
Jo 



dx. 
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and thus the KSMJ metric comes from the metric on A given by 
(53) «u,u>X = J (u'{s),u'{s))ds. 

Again we note that KSMJ were interested in periodic curves on rather than 
odd curves on [—1,1], but the formulas are generahy quite similar apart from 
normalizations. The metric ((53)) has also been studied by Younes et al. |YMSMj : 
a similar metric arises in the study of the Hunter-Saxton equation as well (see 
Khesin-Misiolek |KM| and Lenells [fce]). 

In a sense then, our metric (jS]) lies between the metric ([53| (for which there are 
unique minimizing geodesies and a nondegenerate distance) and the metric (for 
which geodesies cannot be minimizing and the distance is always degenerate) . Our 
geodesies fail to be minimizing by Theorem 14. 2 [ but our distance is nondegenerate 
by Proposition 14.31 Furthermore the geometry induced can be approximated by 
finite-dimensional objects, as in |Pr| . where unit-speed curves are well-approximated 
by a chain of points joined by rigid rods of fixed length, which may be helpful for 
numerical approximations of curves in this geometry. 

Appendix A. Other boundary conditions 

In this paper we have exclusively studied the boundary condition corresponding 
to a whip with one fixed end and one free end. This is the most physically relevant 
condition for an actual whip (a person swings the whip to give it an initial position 
and velocity, then holds the handle basically fixed while the other end swings freely) . 
As shown in |Pr| , the easiest way to handle the technical complications of a fixed end 
is to extend the curve through the origin to be odd; then the boundary conditions 
work out automatically. Hence we have essentially reduced the situation with one 
fixed and one free end to the situation with two free ends. There is no substantial 
difference in any of the results when dealing with two free ends even if the curve is 
not odd. With two fixed ends, the situation is more complicated. (Physically this 
might represent a jump rope being held at both ends.) The same technical issues 
arise, but now it is less obvious how to extend the whip to be odd on both ends; of 
course it can be done, but then we end up with an infinite string and lose some of 
the benefits of compactness. 

Geometrically, however, none of these boundary conditions arc nearly as relevant 
as the periodic condition, since we are interested in curves that form boundaries of 
planar objects (and hence cannot themselves have a boundary). Many aspects of 
this situation are technically easier than the one- fixed/one- free condition we have 
considered, since we can do everything in terms of ordinary Sobolev spaces on the 
circle rather than weighted Sobolev spaces on the interval. The major differences 
are in the upper and lower bounds of the tension, and in the fact that periodicity 
forces 7'(s) ds = 0, which shows up as an extra constraint in some equations. 

Throughout this appendix we work with the circle of length 1. The space of 
curves is the ordinary Sobolev space /C™(S'^) = i7™(S'^,M^), and the subset of arc- 
length parametrized curves is A"HS^) = {76 /C™ : \-f'{s)\^ = 1}. Let T"'-'^{S^) 
denote the space of real- valued functions of class iJ™^^. The proof that A"''{S^) 
is a submanifold of /C™(S'^) is both simpler and works in more cases for periodic 
boundary conditions than for one fixed and one free end. 

Theorem A.l. Ifm ^ 2, then A^'iS^) is a smooth Hilbert submanifold oJJC^{S^). 
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Proof. Define J: JC'" T""-^ by the formula J{ri){s) = \ri'{s)\^. Since H™"! 
functions arc closed under multiplication for rn ^ 2, it is easy to sec that J is well- 
defined and . The differential is easy to compute: we have {DJ)r){u) = 2(u', rj'}. 
Hence for any 7 e A"^{S^) we can prove it is surjcctive; let / e jrm-i 
real-valued function. Write j'{s) = (j){s)dx + '4'{s)dy for functions cj) and V' 

satisfying cfP +1/:'^ = 1, and define u so that u' = {fcp — gijj) dx + {f'ip + 94') Sy, where 
g = acj) + bij: and a and h are constants chosen so that w'(s) ds = (which is 
necessary for u to be periodic). It is easy to see that such constants can always be 
chosen as long as (j) and are not constant, and the only way that could happen is if 
7 were a geodesic. But there are no closed geodesies in M^o Clearly u constructed 
this way is in and so DJ is surjective at any 7 e A^"'{S^). Hence if is the 
constant fimction 1 on 5^, then = A"^[S'^) is a Hilbert submanifold. □ 

The tangent space T^A'"' still consists of 7?™ vector fields u with (u', 7') = 0, so 
that the (formal) orthogonal complement is still 



(T^^™)^ = |A(a7'):a6^"+i| 



(Note that as in Proposition 13. 2[ this only makes sense in T-^/C™ if 7 e A"^'^'^.) So 
the orthogonal projection is P-yiz) = z — [aj')', where 

(54) a" - \j'fa = (z' , 7'), a(0) = a(l), a'(0) = a'(l). 



To prove equation ([54]) always has a solution, we construct the Green function for 
it. 

Proposition A. 2. Suppose 7 and z are smooth. If \y\ is not identically zero, then 
equation (|54p has a unique solution a . 

Proof. Let G: ^ denote the Green function, satisfying 

Gss{s,x) - \i\s)fG{s,x) = -6is-x), G(0,x) = G{l,x),Gs{0,x) = G,(l,x). 
If we can find the Green function, then the solution of (1541) is given by cr(,s) = 

-lo G{s,x)(z'{x),Y{x)}dx. 

Translating by x, we easily see that G{s,x) = ip{s) where 

(55) ^"{s) = K{sf^{s), ^(0) = ^(1), ^'(0) + l = ^'(l), 

with k(s) = \y{s + x)\. So we just need to prove that the boundary value problem 
(|55p has a solution if k is not identically zero. 

Let ipi and denote the solutions of (j55p with boundary conditions V'i(O) = 1, 
ip[{0) = 0, V32(0) = 0, and ^^^(0) = 1. Clearly ipi{s) ^ 1 and ^'2(3) ^ 1 for aU s. 
We can write (p = Aipi + Bif2 where A and B satisfy 

(56) A[cri(l) - 1] 5(72(1) = and AcrJ (1) B[4(l) - 1] = 1. 

Using the reduction of order trick, we can compute that (72 (s) = o'i(s) % dr/ai(ry, 
so that the determinant of the system ([55]) is 



1 [a,{l)-lf-a[{l)Cdx/a,{xf, 
Jo 



■^The same construction works for the arc space on any Riemannian manifold M, and A"^(S^) 
{7: 5^ — > M : |7'| = 1} fails to be a manifold if M contains any closed geodesic of length 1. 
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which is always negative unless cri(l) = 1 and <7'i{l) = (which happens if and only 
if K = 0). Hence we can solve for A and B, and so we obtain the Green function. 
Uniqueness is trivial using a standard energy argument. □ 

Now we obtain upper and lower bounds for a. 
Proposition A. 3. The solution ip of (|55p satisfies 

(57) e"^/^/e ^ V9(s) ^ 1 + , where g = ^ K^{x)dx, 

for all s e . 

Proof. To prove the lower bound, set 7 = hiip, so that 7" = — 7'^. Let ti be a 
point where 7' = 0; since 7 is convex, "f{u) is the minimum of 7. Integrating by 
parts twice from to m we get j{u) ^ 7(0) — xk'^{x) dx, and similarly integrating 

from u to 1 we get 7('ii) ^ 7(1) — Jg(l — x)k^{x) dx. Since 7(0) = 7(1), averaging 
these estimates gives 



(0) - 1 f.^ 
Jo 



7(u) ^ 7(0) - i I K\x)dx. 

Thus we get ip(u) ^ ip(Q)e~Bl'^. 

Now we need a lower bound for <p(0). Since </?'(!) = <y3'(0) + 1 and = 
we have 7'(1) = 7'(0) + 1/(^5(0), so that 

= 7'(1) - 7'(0) = £ n\s) ds - £ j'isf ds ^ g. 

Combining this with our estimate for ip{u) gives the lower bound (p{u) ^ e~^/^/p. 

The upper bound is obtained differently. Writing if = Atpi + Bip2 as in Propo- 
sition |X2j we get 



L^'(l) + [v.i(l)-lM(l)]2' 

where L = dx / ipi{xY . Hence A ^ ipi{l) / ip'i{l) . Now since Lpi is convex, we know 
(^']^(1) ^ <(2i(l) — 1, which implies ^ ^ 1 + l/(y9'(l). Finally, we use the fact that 
fi{s) ^ 1 for all s to get 

<^;(1)= J Lp-^{s)K^{s)ds^ ^ K{sfds^ \j"{s+x)\ds^ > ^TT^ , 

using the well-known bound on the total curvature of a closed curve in the plane 
(see for example do Carmo |dCl| . Section 5-7). □ 

As a corollary, we get a bound for the curvature in exactly the same way as 
Theorem 15.11 

Proposition A. 4. The sectional curvature of A{S^) in the metric is bounded 
below at J e A{S^) by 47r^e^^/^/p, where g = k(s)^ ds, and is unbounded above. 

Proof. The proof is the same as Theorem 15.11 The factor 47r^ comes from the fact 
that the smallest eigenvalue of the derivative operator has size 2n. □ 
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We can prove using techniques similar to |Pr| the analogue of Theorem 13.41 
that solutions of the geodesic equation exist as long as 77(0) 6 A^{S^) and 77(0) 6 
T~^A'^{S^); the weakening comes from the fact that we never need to use the esti- 
mates of Lemma |2. II to fix the weighting at the endpoints. 

Unboundedness of the curvature again implies that the exponential map cannot 
be C^; this can be checked explicitly using simple explicit solutions. The simplest 
explicit solution of the geodesic equation on is 

(58) Tj{t,s) =-f{s+ujt), uj{t,s)=uj\ 

where 7: S*^ ^ is any closed curve. We can compute Jacobi fields along such 
a curve explicitly (for example, if 7 is a circle) and show that they have zeroes for 
arbitrarily short times, as in Theorem 14.21 

One might object to the notion that unphysical solutions such as ([58)) should 
be allowed, especially in application to shape recognition, since the image of the 
curve doesn't change with time. The typical way to resolve this (see for example 
[YMSMj ) is to quotient out by the translations. The metric on -4(5'^) is obviously 
invariant under the action by translations, so we get a metric on the quotient 
A{S^)/S^. Since ([58)) is always a geodesic, this is a Riemannian submersion with 
totally geodesic fibers. This corresponds to requiring tangent vectors to satisfy 
not only (u',7') = but also ^gi(u{s),Y{s)} ds = 0. The new orthogonal space 
Tj{A{S'^)/S^) then consists of fields of the form -j^icrj') + cj' for functions a and 
constants c. The geodesic equation is then 

rjtt = dsicrrjs) +ajs, cTss - {ilssl'^o- = -\ristf, 

and we have ^ \gi(j]t,'nsyds = c, so that c = in order to preserve horizontality. 
Hence the same geodesic equation guarantees that the translations disappear as 
long as ^gi(r]t ,r]s} ds = initially. 

Appendix B. Removing the length constraint 

We note that one drawback to these equations in shape analysis is that the space 
of all odd curves (modulo reparametrizations) is not exactly the same as the arc- 
length parametrized curves, since the space of all curves includes those of arbitrary 
length while ours consists only of curves of length 2. To extend this, we would have 
to work with a slightly different version of A. Let 

:4 = {7 : [-1, 1] ^ r2 : ly I = const and 7(-s) = -7(3) Vs} . 
Theorem B.l. The geodesic equation on A in the weak metric ([8]) is given by 

(59) vtt = ds{<J7j,), \Vs\^=e^ 

(60) fass-\ilss\'' + \Vst\''=C, J <jds = 0, a(-l) = a(f) =0. 

Equation (|60p always has a solution for any given rj. Here £ and C are constant in 
space but not necessarily in time, and we have ^{P) = 2C. 

Proof. The condition for a vector field u along 7 e ^ to be in T^A is that j^(u', 7') = 
0, so that the orthogonal complement consists of fields of the form ^(177') for func- 
tions a with (t(s) ds = 0. The geodesic equation therefore still has the form of 
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the wave equation ([59| . Here \r]s\ = I for some I which is constant in space but 
may depend on time. 

To find the equation for a we differentiate the constraint i5s|?7sp = twice with 
respect to time to get d a((j\stt^r\s^ + \r\at^) = 0. Plugging in rju from ([59]) . we obtain 

-|77s.f + =0, 



which integrates to (|60| . 

To prove we can actually solve (|60p . we consider the slightly modified Green 
function Gi(s,x) satisfying 

fd^,G{s,x) - \vss{s)f = -S{s - x), G(-1,.t) = G(l,x) = 0. 

If we write the solution a in terms of this Green function, then C is determined by 
the fact that the integral of a vanishes, i.e., that 



J J G{s,x){\vtAx)\'' -C)dxds = 0. 



By Theorem 13. 11 G(s, x) ^ for all s and x in the square [—1, 1]^ and is zero only 
on the boundary, which means we can always solve this equation for C. It is easy 
to compute that ■£^(e'^) = 2C. □ 

Although we can solve (|60p , it is far from clear that we can solve (|59p : the diffi- 
culty is that since a integrates to zero, it cannot be strictly positive, which means 
([5^ is always of mixed type, and hence substantially more difficult to analyze. We 
also lose positivity of the sectional curvature as in Theorem 15.11 In addition, the 
fact that C is always positive except in degenerate cases means that i is always 
increasing, so that in particular A is not a totally geodesic submanifold of A. We 
leave these complicated issues aside for now however, noting merely that the situ- 
ation is not at all improved by changing the boundary conditions as in Appendix 

m 
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